Polaron relaxation in a quantum dot due to anharmonic coupling within a mean-field 

approach 
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We study the electronic relaxation in a quantum dot within the polaron approach, by focusing on 
the reversible anharmonic decay of longitudinal optical (LO) phonons forming the polaron into lon- 
gitudinal acoustic (LA) phonons. The coherent coupling between the LO and LA phonons is treated 
within a mean-field approach. We derive a temperature-dependent inter-level coupling parameter, 
related to the Griineisen parameter and the thermal expansion coefficient, that characterizes an 
effective decay channel for the electronic (or excitonic) states. Within this theory, we obtain a 
characteristic anharmonic decay time of Ins, 2-3 orders of magnitude longer than previous predic- 
tions based on the Fermi's Golden Rule. We suggest that coherent relaxation due to carrier-carrier 
interaction is an efficient alternative to the (too slow) polaron decay. 

PACS numbers: 71.38.+I; 73.61. Ey 



CD 



O 

o 



> 
o 
in 

o 

O 



X 



Introduction. The understanding of the relaxation 
mechanisms in quantum dots (QDs) is important from 
both, the fundamental and technological points of view 
since QDs are increasingly finding applications in lasers, 
single-photon emitters, and possibly quantum computers 
where fast carrier dynamics are indispensable. Whereas 
in bulk semiconductors, the hot carrier relaxation is me- 
diated via emission of optical phonons^, in a QD the dis- 
crete electron and hole levels are separated by energies 
that generally do not match any phonon energy. This 
simple idea gave rise to the " phonon bottleneck" concept 
which predicts the inefficiency of hot carrier relaxation 
by emission of phonons in QDs^. However, this predic- 
tion relies on the assumption that the phonon emission is 
irreversible, with a probability described by the Fermi's 
Golden Rule (PGR). 

It has then been realized that multiple scattering 
processes are important and that the electron-phonon 
(e-ph) interaction in QDs must be treated in a non- 
perturbativeS. and non-adiabatic"* way, leading to the en- 
ergy spectra described by polaronic quasi-particle excita- 
tions. Polaron spectra of QDs have been considered theo- 
retically using several approaches^'^'^ and various steady- 
state observable properties [such as photoluminescence 
(PL) and Raman spectra] have been calculated and com- 
pared to experimental data, convincingly demonstrating 
the importance of the polaron concept^'^'^ i ^1^1^° . 

One might then question whether the "phonon bot- 
tleneck" in QDs really exists and, in fact, several experi- 
ments point in either direction. On one hand, an efficient 
relaxation of optically created electron-hole pairs (here- 
after called excitons) was reported in a number of works 
studying self-assembled QDs (SAQDs)iiii^ii^, with both 
PL rise tim e ^-'^'l" and photoinduced intraband absorp- 
tion decay time^^ below or of the order of 10 ps. Recent 
studiesi^ii^ performed on chemically grown nanocrystals 
(NCs), where exciton energy level spacings are larger 
than in SAQDs, also revealed ultrafast intraband relax- 
ation, although the mechanism seems to be different for 
CdSe and PbSe NCs. There are some works where the 



relaxation of a lone carrier (e.g. electron in an n-doped 
QD) was studied^iiiiii^. Their results suggest that the 
relaxation is slower than for excitons but still fast, with 
a characteristic time of several tens of picoseconds for n- 
doped InAs/GaAs SAQDs^iil. On the other hand, there 
are published experimental results that support the ex- 
istence of a phonon bottleneck effect in the relaxation of 
optically created excitons, both in self-assemble d^^'^° and 
nanocrystal^l'j^ QDs. For instance, a relaxation time of 
7.7 ns, ~ 15 times the radiation lifetime, was obtained 
for InAs/GaAs SAQDs in Ref.[li. 

The polaron model can explain intraband relaxation 
of carriers in QDs only within the context of so-called 
pseudo-relaxation^- (i.e. oscillatory dynamics), since po- 
larons are stationary states of an electron (exciton) cou- 
pled to optical phonons. Some additional interactions 
should therefore be responsible for the true polaron relax- 
ation (i.e. thermalization) . Several possible mechanisms 
of hot carrier relaxation in QDs have been proposed: 

(i) The polaron has a rather short lifetimeSi because of 
the anharmonic effects that lead to a fast decay of con- 
fined optical phonons forming the polaron. This mech- 
anism has been considered for both excitor^^ and lone 
electron^Si^l relaxation. 

(ii) In SAQDs, the polaron can relax via an Auger-type 
mechanism assisted by electrons present in the wetting 
layer where the energy spectrum is continuous^*'^^. 

(iii) Acoustic phonons can provide the possibility of 
transitions between different (exciton-) polaron states^. 
If the acoustic phonon spectrum is (at least partially) 
continuous and the polaron spectrum is sufficiently dense, 
this interaction would drive the polaron dynamics to- 
wards equilibrium. 

An Auger-type mechanism was also proposed'^*'''^^ for 
the relaxation in chemically grown QDs, according to 
which the excess energy is first transferred from the elec- 
tron to the QD hole through their Coulomb interaction 
and the subsequent hole cooling occurs via emission of 
acoustic phonons because the hole level spacings are rela- 
tively small and match the continuum of acoustic phonon 
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energies. Experimental data support this mechanism, 
at least for CdSe nanocrystals,i4i^ii^ If one uses the 
exciton-polaron language, this mechanism is equivalent 
to (iii). 

In this Brief Report, we show that the anharmonic- 
ity mechanism (i) is too slow and therefore not rele- 
vant in most of the experimental situations. It was ini- 
tially proposed and illustrated for the simple case of 
a single longitudinal optical (LO) phonon mode cou- 
pled to two electron states in Ref. [2J using the rotating 
wave approximation. Later, more realistic calculations 
of the polaron spectrum and relaxation rate were per- 
formed for an exciton— and a lone electron^SiSl in an 
InAs/GaAs SAQD. In the latter worls^l, a good agree- 
ment with experimental findings of Ref. 1/ has been 
achieved by including several anharmonic phonon decay 
channels known in the literature for various materials^^. 
However, in all these works the LO phonon decay was 
considered to be irreversible and FGR was used to cal- 
culate the polaron relaxation rate. We shall avoid this 
approximation, which is known to be a rather crude one 
for QDs and incoherent with the polaron concept. In- 
stead, we develop a non-perturbative mean-field theory of 
the QD polaron coupled to acoustic phonons through re- 
versible anharmonic processes, resulting in a much slower 
relaxation. Our main conclusion is that the mechanism 
(i) is irrelevant in most experimental situations. 

The model and approach. We present a non- 
perturbative approach to the non-equilibrium problem 
of electron relaxation in a quantum dot. This shall be 
discussed within a simple (minimal) model which nev- 
ertheless includes all relevant interaction terms. Our 
QD contains two non-degenerate electronic levels (sep- 
arated by an energy e), coupled to a single confined LO 
phonon mode. The LO phonons can reversibly decay into 
a couple of longitudinal acoustic (LA) phonons through 
an anharmonic process that has been identified experi- 
mentally for some bulk semiconductors'^^. We shall also 
include the possibility of further anharmonic interactions 
of these (secondary) LA phonons with other vibrational 
modes in the QD or its surroundings. The Hamiltonian 
of this model can be written as follows: 

H — Hpol + Hla + Hint + Hanharm (1) 

where the polaron Hamiltonian is given by 

Hp = €c\ci+ujBb'^b + M{c\ca + clci){b + b'^) , (2) 

the Hamiltonian of the acoustic phonons reads 

Hla = ujAO^a , (3) 

and the interaction term is 

Hint = G{a^a% + b^aa) . (4) 

In Eqs. ujb and lua are the energies of the LO and 

LA phonons, respectively (we set h = 1), and M is the 



inter- level e-ph coupling constant^. In Eq. ([4]), G de- 
notes the characteristic energy of the LO ^ 2LA process 
related to the experimentally measured anharmonic de- 
cay time constant'^^-'^^. The last term in Eq. ^ stands 
for anharmonic interactions of the acoustic phonons (ex- 
cluding those written explicitly as Hint)- This term is 
needed for a well-defined ground state (see below) . 

In the spirit of a mean-field approach, the interaction 
term shall be approximated as 



Hi, 



aa) + a'' a' {b) + (b^}aa) , (5) 



thus neglecting the coupling of the fluctuations of the 
bosonic fields. By defining shifted bosonic modes, b 
b + G{a^a^) /(jJb, the Hamiltonian ([1]) can be written as 



(G(atat))2 
lub 



H — Hpol 
where we have 
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LA 
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-M{c\co + cjci) 



and 



Hla — ujao-^o- + D{a^a^ 



(6) 



(7) 



(8) 



In Eq. (O , we have introduced a new inter-level coupling 
constant. 



2G'(atat^ 
lub 



M = M 



(9) 



and in Eq. ^ an effective coupling constant D including 
all the anharmonicity effects involving the LA phonons. 
This constant takes into account the effect of the last two 
terms in Eq. ([5]) as well as Hanharm which represents 
fourth-order processes with the participation of two LA 
phonons, i.e., D = G{b) + D. 

There is a new direct relaxation term A_ff , character- 
izing the new relaxation scale r = 1/M. To quantify this 
scale, we need to know the value of the anomalous am- 
plitude (a^flt). The explicit evaluation of the anomalous 
amplitude involves the diagonalization of Hla, which can 
be made using a Bogoljubov-Valatin transformation^, 
1] = ga + ha\ with the normalization of the canonical 
commutator relations ~ — 1. This gives 



g = —T 1 h = 

USA n 



i^A - i^A 
2UJA 



with the renormalized eigenfrequency 
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The occupation number of the acoustic phonons can 
be expressed in terms of the expectation value of the 



new bosons, ^r\^r\) 
riA = (a^a) = ribc 
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a^a) = ribos + h'^{l + 2nbos) 
amplitude entering Eq. ([9]) is given by 

(aa) = (a^a^) = ~gh{l + 2nb. 



The anomalous 
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FIG. 1: (Color online) EfTective anharmonic inter-level cou- 
pling parameter Af calculated for excitons and lone electrons 
in CdSe and InP nanocrystal QDs of 2 nm in radius. 



Thus, all interaction terms in Eq. ^ (except Hpoi) are 
now expressed in terms of the effective coupling constant 
D. In particular, the new inter-level coupling constant 
^ is given by 



2GD 

M = -M- {l + 2nbos) 



Self- consistent mean field theory. In a self-consistent 
mean-field theory involving only electrons and LO and 
LA phonons, we should put D = G{b). The self- 
consistency condition is then defined by the following two 
equations: 



{b) = 
{aa) = 



2G( 



lub 



-G{b) 



(14) 

-,{l + 2nbos). (15) 



This yields the following self-consistent equation for the 
anomalous amplitude: 

Luluj% - WC^iaaf = AG^{1 + 2nbosf ■ (16) 

For G — > 0, we thus have: 

(aa) « c^awb/(4G2); Mmf ~ Mw^/(2G). 

This means that the characteristic decay time, t ^ 
1/Mmf, tends to zero for vanishing anharmonic cou- 
pling, which certainly is an artifact of the model. The 
physical reason for this is the omission of other anhar- 
monic processes [D — 0). It is well known that self- 
consistent mean-field theories for boson-fermion systems 
may lead to such inconsistencies!^ 



Effective anharmonic coupling constant. In order to 
overcome the difficulty of the previous section, we have 
to take into account the fact that LO 2LA is not the 
only anharmonic process involving the LA phonons con- 
sidered in our model. In Eq. ([8]) we have introduced 
D as an effective coupling constant representing all the 
third- and fourth-order processes with the participation 
of two LA phonons. It will now be considered indepen- 
dent of G, which is the coupling constant of one partic- 
ular anharmonic process. Assuming that we can ignore 
other cubic anharmonicity effects (involving just one LA 
phonon), we can relate D to the experimentally known 
Griineisen parameter—, — ^dhiO/dhiV where is 
the Debye temperature and V the crystal volume. In the 
Griineisen approximation, the variation of all vibrational 
frequencies with the crystal volume is assumed to be the 
same^, dh\0 ~ d\nuj. If we consider the effective cou- 
pling constant as a function of temperature, D(T), and 
use Eq. (jlip for uja, then we can write: 



din ua 



1 dlnujA 



oiT dT 



1 



2 dD^ 
Q\~dT 



where ax is the thermal expansion coefficient. Thus, we 
obtain the following equation for D{T), 

2 dD^ 



iGrOLT 



(13) Its solution, with the condition -D(O) = 0, yields: 



2 ^ 



exp {~2-iQraTT)\ 



1/2 



(17) 



(18) 



With typical values of olt ~ 10 "^K ^ and ^gt ~ 1-5, 
we obtain Dj^A ~ 1/20 for T « 300K. Accordingly, 
ujA ~ L^A- Fig. [1] shows the temperature dependence of 
the effective inter-level coupling parameter M following 
from Eqs. and ([T5)) . calculated for spherical QDs. 

We assumed wa — '-^b/'2 and used the anharmonic de- 
cay constants from experiment (G=0.15 ps^^ for CdSe— 
and 0.04 ps~^ for InP—) and calculated e-ph couphng 
parametera^i^S. M thus turns out to be very small (even 
for electron in CdSe QD), much smaller then the typical 
e-ph interaction energy (w 3 meV for excitons and ~ 20 
meV for electrons interacting with spherical-symmetric 
phonons). It is larger for lone electron, compared to exci- 
ton, because of the stronger coupling to polar phononsj^ 
To conclude, even though the anharmonic decay of LO 
phonons is faster in QDs than in bulk semiconductors 
because of the absence of the restrictions due to momen- 
tum conservation, our results indicate that the energy 
scale of this interaction is far below the other relevant 
characteristic energies. Consequently, this effective de- 
cay channel is too slow to account for the experimental 
data. On the contrary, coherent coupling between carri- 
ers can produce effective (pseudo-) relaxation of the QD 
polaron, with r ~ 10 ps, even if there are no dissipative 
relaxation channels present. 
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Summary. Assuming the coherent couphng between 
the optical and acoustic phonons and using a mean-field 
approximation, the Hamiltonian ([1]) describing the inter- 
acting polaron can be reduced to H ~ Hpoi -\- AH with 
AH given by Eqs. 0, IHl) and Note that AH 

is temperature dependent and [apart from the parame- 
ters of the "non-interacting" Hamiltonian entirely 
determined by the thermal expansion coefhcient and the 
Griineisen parameter. The term AH corresponds to a 
(pseudo-) relaxation time r ~ 1 ns, i.e., 2-3 orders of 
magnitude larger than the previous predictions based on 
the Fermi's Golden rule. This renders the anharmonic 
decay channel virtually irrelevant, at least for nanocrys- 
tal QDs where the polaron is formed by a small number of 
confined optical phonon modes most strongly coupled to 
the electron or exciton. Instead, the temporal evolution 
of the photoinduced polarization in a QD can proceed 
much faster just because of the quantum beats result- 
ing from coherent many-body interactions ,22. In doped 



SAQD heterostructures, the relaxation of the lone elec- 
tron or hole occurs due to the coupling to the carriers 
in the wetting layer,— while in optically excited QDs the 
same mechanism operates because of the electron-hole 
interaction^^. In some cases, this pseudo-relaxation can 
be taken over by the true (dissipative) one mediated by 
the continuum of acoustic phonons.-'' In fact, atomistic 
calculations using the time-domain ab initio approach in- 
dicate that the relaxation in nanocrystals occurs via e-ph 
coupling to acoustic-type vibrations and the phonon bot- 
tleneck can exist only at low energies in very small dots.— 
We finally note that recent experiments^ have shown a 
bottleneck effect on the hole relaxation in CdSe QDs of 
1-3 nm in radius (with the gap between the ground and 
first excited states w 30 - 100 meV, similar to SAQDs), 
which is consistent with the conclusions of this work. 
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